Theory of Coherent Raman Superradiance Imaging of Condensed Bose Gases 
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We model the ofF-resonant superradiant Raman scattering of light from a cigar-shaped atomic 
Bose-Einstein condensate. Absorption images of transmitted light serve as a direct probe of long 
range coherence in the condensate. Our multimode theory is in good agreement with the time- 
dependent spatial features observed in recent experiments, and the inclusion of quantum fluctuations 
in the initial stages of the superradiant emission accounts well for shot-to-shot fluctuations. 



I. INTRODUCTION 

In recent work, L.E. Sadler et al. discuss the use of 
superradiant Raman scattering as a probe of long range 
coherence in an ultra-cold bosonic gas 0. Their inves- 
tigation is one of few experimental studies of extended 
sample superradiance that is not only temporally, but 
also spatially resolved. 

Superradiance is a well-known phenomenon |2| that 
was first discussed by Dicke Q in 1954. In that process, 
atoms scatter light collectively, producing a short burst 
of intense radiation. In a typical superradiance experi- 
ment with Bose-Einstein condensates, an elongated sam- 
ple of ultracold atoms is subjected to a pump pulse of 
far off-resonant laser light. For an appropriate polariza- 
tion of that field, the cigar-shaped condensate geometry 
results in the scattering of photons predominantly into 
modes propagating along the long axis of the sample, 
commonly referred to as the end-fire- modes (EFM). As 
the atomic cloud is Bose-condensed, the scattering pro- 
cess simultaneously leads to the coherent amplification 
of the recoiling atomic fields. Such coherent matter- wave 
amplification (CMA) has been observed both in the case 
of Rayleigh scattering ^ ^, |6j| and of Raman scattering 
0, , and several theoretical descriptions of this effect 
have been published "o!, [HITl [H . 

A necessary condition for superradiance to occur is 0, 

HI HI 



(1) 



where T2 is the reciprocal inhomogeneous linewidth and 
Tc = 2/(cp7A^)^/^ is the superradiant cooperation time 
with p the atomic density, c the speed of light, A the 
wavelength and 7 the linewidth of an isolated atom. 

An ultracold Bosonic gas at temperature T with < 
T < Tc, where Tc is the critical temperature, consists of 
both condensed and non-condensed phases. As Doppler 
broadening for the two phases can be dramatically dif- 
ferent, it is possible for the appropriate choice of experi- 
mental parameters to have condition Eq. ^ satisfied for 
the condensed phase, but not the non-condensed phase 
0. Under such circumstances light will scatter super- 
radiantly only from the condensed portion of the gas, 
thus providing the experimenter with a sensitive probe 
of that phase. In Ref. [H this approach is followed to 



both quantify the condensate number and study the spa- 
tial and temporal evolution of the superradiant process. 
It provides a novel probe of long-range coherence which 
may lead to insight in the symmetry-breaking dynamics 
of normal- to superfluid phase transitions. To our knowl- 
edge only one previous study [3| exists that resolves spa- 
tial features in extended sample superradiance in atomic 
vapors. In that case the angular pattern of the emitted 
EFM was integrated over the course of a superradiant 
pulse and found to consist of several bright spots (Fig. 
3(a) in Ref. 4]). 

Zobay and Nikolopoulos [T^ [l7j have given a detailed 
semiclassical analysis of the spatial features of both the 
matter-wave and the optical fields in CMA experiments 
based on Rayleigh scattering. They found that propaga- 
tion effects play a crucial role in the amplification process 
and account for several characteristic features seen in ex- 
periments 0, Q- These include the characteristic "X" 
and fan shapes of the atomic recoil modes correspond- 
ing to the strong and weak pulse limits respectively, the 
asymmetry between forward and backward side-modes in 
the strong pulse regime, and the depletion of the conden- 
sate center in the weak pulse regime. Furthermore, they 
found that low superradiant emission does not necessarily 
imply a small EFM field inside the atomic sample, as a re- 
sult of the scattering of EFM photons back into the pump 
beam. They also demonstrated that propagation effects 
lead to sub-exponential growth of the scattered fields in 
contrast to the exponential growth seen in fully quan- 
tized uniform field models, an effect reminiscent of laser 
lethargy in short- wavelength optical amplifiers 0, 0| . 

An important aspect of superradiance experiments is 
the appearance of large shot-to-shot fluctuations, a re- 
sult of the quantum noise that dominates the dynamics 
during the initial stages of the experiment [2, 120I |2]1 . 
At later times the scattered fields become macroscopi- 
cally occupied and evolve in an essentially classical man- 
ner from initial conditions determined by those quantum 
fiuctuations. A common strategy to treat the full evolu- 
tion of superradiance is therefore to break the problem 
into an initial quantum stage followed by a classical stage. 
The initial stage is analyzed primarily to determine the 
appropriate probability distribution of initial conditions 
for the classical stage. Experimentally observed fluctua- 
tions in the classical stage are simulated by solving the 
dynamics a large number of times with initial conditions 
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FIG. 1: Three-level Raman transition 
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FIG. 2: Experimental setup. 



chosen according to this probabiUty distribution, each 
simulation representing a single realization of the exper- 
iment. 

The primary aim of this paper is to provide a de- 
tailed model to complement the experiments reported in 
Ref. [3. We generalize the semiclassical model of Refs. 
0, llTj l to the Raman scattering case, and extend it to 
a multimode quantum description that allows for a sys- 
tematic treatment of the build-up of the classical fields 
from quantum noise, thereby accounting for shot-to-shot 
fluctuations. Whereas Refs. 0, ^3 focussed on the 
post-pump expansion patterns of the recoiling modes to 
check self-consistency of their predictions, we investigate 
the time-dependent imaging of the condensate while un- 
dergoing superradiant emission. We find good qualitative 
agreement with experimental observations. 

The paper is organized as follows. Section ^ dis- 
cusses the experiment under consideration, introduces 
our model, and considers the initial stages of the evo- 
lution of the atoms and light field, treating the scattered 
optical field quantum mechanically. Section IIIII turns 
then to the classical stage by solving coupled Maxwell- 



Schrodinger equations within the slowly varying envelope 
approximation. Numerical results are presented in Sec. 
ITvl and Sec. E\ is a summary and conclusion. 



II. INITIAL STAGE AND QUANTUM 
FLUCTUATIONS 

We consider a cigar-shaped Bose-Einstein condensate 
of width w and length L with a total of N atoms. The 
ultra-cold atoms in the condensate undergo Raman scat- 
tering between two ground states |1) and |2) via an off- 
resonant excited state |e). The energy of level |1) is taken 
as the zero of energy, and the energies of the states |2) 
and |e) are hui2 and hcoe, respectively, see Fig. 1. We as- 
sume that the transition |1) |e) is driven by a classical 
pump laser Ei(t) of frequency lul and polarized along 
the y-axis, see Fig. 2, with 



(2) 



while the transition |e) — s- |2) takes place via spontaneous 
emission to a continuum of vacuum modes with frequen- 
cies 



Wfc = We - W2 - 4- 

The total electric field is then 



(3) 
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h.c. 



where the creation operators aly^{t) obey the usual 
bosonic commutation relations 



= (5(k-k')fe'. 



(4) 



The incident laser field envelope EL{t) is taken as con- 
stant in amplitude during the initial stages of the am- 
plification process, but its full time dependence can be 
accounted for in the classical stages of the evolution. 

We proceed by introducing bosonic matter-field cre- 
ation and annihilation operators '0|(r,f), ('0i(r, i)), that 
create (annihilate) an atom at position r in electronic 
state \i) = |e) or |2), with 



Mr,t),^]{r',t) = <5,,(5(r - r') 



(5) 



in terms of which the Hamiltonian of the atom-field sys- 
tem is H = Hq + He, with 
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(6) 



and the interaction Haniiltonian 



Hr 



(7) 



describes the electric dipole interaction between the 
atoms and the electromagnetic field, being the dipole 
moment of the \i) |e) transition. 

We assume that the atoms are initially in their ground 
state |1), and that the pump laser is sufficiently far de- 
tuned from resonance for the excited state population to 
remain negligible at all times. The excited level |e) can 
then be adiabatically eliminated in the standard fashion. 
Thus transforming to slowly varying interaction picture 
operators 



(8) 



^i(r,t) = ?/'i(r,t) 



(9) 
(10) 

(11) 



and performing the rotating wave approximation (RWA) , 
yields the effective Hamiltonian 



dr / dk 



{.9(k)atj^^^(r,i)e'(''^-'')-^^i(r,t)e^(*^'*^)* + /I.C.} , 



(12) 



where 
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(13) 



In Eq. (|12(l we have neglected higher-order scattering 
terms coupling modes of the scattered light field, an ap- 
proximation justified in the early stages. 

The effect of the Hamiltonian H is to transfer 
atoms initially in the ground state |1) and with spa- 
tial wave function 0o(r) into the ground state |2) 
with a momentum-conserving spatial wave function 
(j)o{r)e^''^''~^'>''^ via scattering of a photon into mode k. 
Following Ref. 10] we therefore expand the matter- field 
operators into quasi-modes according to 



V'i(r) = 0o(r)co 



and 



Mr)^ / dq<^o(r)e''i" 



(14) 



(15) 



where cq annihilates a particle in a quasi-mode with elec- 
tronic state |1) and wave function (r|0) — 0o(r) and 
Cq annihilates a particle in a quasi-mode with electronic 
state |2) and wave function (r|q) = (j)o{r)e^'^''^ . These 
quasi-modes are nearly orthogonal provided that they are 
separated by an angle 

> 2K^/\k\ = X/inw), 

where = 2/w is the momentum width of the conden- 
sate. In that case the quasi-mode creation and annihila- 
tion operators Cg obey to a good approximation bosonic 
commutation relations 



Sqq' 



(16) 



and are statistically independent during the early stages 
of the evolution. 

In terms of the quasi-modes q, the Hamiltonian 1)12(1 
becomes 

Hc^Y.j'^^ {?7(k, q) e^('^-''=)*44^o + h.c^ (17) 
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where 



r?q(k)=5(k) / dvlMrre 



|2 i(kL-k-q)T 



(18) 



We have used the normahzation condition J dr \(j>o{r)\'^ — 
1. The Hamihonian (|17|l differs from the effective Hamil- 
tonian derived by Moore and Meystre for the Rayleigh 
scattering case [I^ in that the electronic state of recoiling 
atomic modes is |2) instead of the original ground state 
|1) and the term oscillating with frequency Sl — Sk = 
LUk — {lul — '-^2) enforces the coupling of the scattered 
modes to photons with frequency lol — ^2 instead of w^. 

With these differences in mind, we may easily gen- 
eralize the results of Ref. [l^ and determine that the 
quasi-modes initially grow exponentially, 

Cq = exp[GqiVV2]£q(0) + f dr exp[Gq7Vt/2]/t(t - r), 

Jo 

(19) 

where 



(20) 



and fq{t — r) is a noise operator the second-order cor- 
relation functions of which are given in the Markov ap- 
proximation by 

{flit)Ut')) = 0, (21) 

imfjit')) = G^NS{t-t'). (22) 

In this limit, the probability Pq{n,t) of having n atoms 
in mode q at time t is that of a chaotic field, 



1 



1 



lit) 



(23) 



where nq{t) = (c^Cq) is the mean number of atoms in 
the quasi-mode q at time t. 

An important feature of the linear gain factor Gq is 
that it remains relatively constant for quasi-modes ex- 
cited via the scattering of photons at small angles 6*1^ 
with respect to the long axis of the condensate. This is 
the case until 0k reaches the geometric angle 

after which the linear gain Gq falls off rapidly. The elec- 
tromagnetic modes corresponding to scattering into that 
angle collectively form the EFMs, and they dominate the 
short-time dynamics of the system. 

This suggests that we may accurately simulate the lin- 
ear dynamics of the superradiant system by considering 
the scattering of photons into a finite number of modes 
distributed within a solid angle 2TT9g around the long axis 
of the condensate only. We note that while the geometric 
angle into which significant scattering takes place is fixed 
by the aspect ratio of the condensate, the number m of 
independent quasi-modes depends also on the wavelength 
of the optical fields involved. 



m 



2\L 



(24) 



where F is the Fresnel number. Typical experiments cor- 
respond to a number of quasi-modes m « 1 ^ 10^. 



III. CLASSICAL EVOLUTION 

We now turn to the quantum-noise initiated classical 
regime that occurs once the supcrradiance process is fully 
underway and the scattered modes are macroscopically 
occupied. In typical experiments of Ref. 0] the velocity 
of the recoiling atoms is such that they may traverse half 
a condensate width over the course of a single run of 
the experiment. In addition, atom-atom scattering may 
cause significant dephasing on the same time scale. In 
this section we therefore include both kinetic energy and 
mean field atom-atom collisional terms. 

In the classical regime the optical field can be described 
by a complex field amplitude 



-i(\iL-T-Ui^t) 



y [i?L(r,i)e'('^-'-"-*)+i?2(r,t)f 



Following the discussion of Sec. ^1 we restrict the sum 
over k to the discrete set of quasi-modes in the end-fire 
cone and set LUk = u>l—'^2- At this point we no longer as- 
sume that the pump field remains undepleted. As usual, 
the field amplitudes Ei{r,t) are assumed to be slowly 
varying. 



V£;,(r,i)| « \hE,{r,t)\, 
« |..^;.(r,t)|. 



dt 



(25) 
(26) 



Eliminating adiabatically the excited electronic state of 
the atoms as before by introducing the slowly varying 
Schrodinger field operators 

Mr,t) = Mr,t) (27) 
V;kjr,i) = ^,(r,t)e^("=+"")*e-'''-- (28) 
M^,t) = V'k(r,t)e'('^^+'^=-^*e-*(''^"'~''-''\ (29) 

where huer = 'h^k\/2m and 'hij02r = ^-^(^i + k'^)/2m, 
results in the effective Hamiltonian H ~ Hq + Haa + He 
where 



Ho 



/-{? 



2m 



(30) 
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(31) 



Hr = 



9ij 



(1 - i-//Sl) 



V'l(r,i)V'j(r,i)- 



(32) 



r 



Here 



(33) 



in which the factor 3 arises from an average over all pos- 
sible orientations of the dipole and 



Uii{kk) 



^U^^^^{a + ika^), (34) 



where a is the s-wave scattering length and the indices i, j 
run over — fc, 0, k. The last term on the right in Eq. 134() 



is a momentum-dependent loss that has been used to 
account for elastic scattering losses in the Gross-Pitaevski 
equation |23| . It arises upon keeping the second-order 
term when expanding the manybody T-matrix in powers 
of k. For the Raman system under consideration we need 
not include higher-order modes of the form V'nkt+qk {n,q 
being integers) in the sum over k, since the second ground 
state can no longer absorp pump photons 0- 

Within the slowly varying envelope approximation 
(SVEA) the Hamiltonian H32|l yields the matter-wave 
Heisenberg equations of motion 



rf^o(r,t) 
dr 

rf^k(r,t) 
dr 



(1 - i^/5l) 



fL*0 + xE^k*-k 

di ^ 



j 



(35) 
(36) 



r 



We have cast these equations in dimensionless units 
where the dimensionless time is 



ol 



in which 



dEf 
V3?i 



(37) 



(38) 



is the the Rabi frequency of the incident laser field of am- 
plitude E^^ and d = ^/dld2 ■ The elecric field amplitudes 
are rescaled as 



E[" 



and the matter-wave fields as 



where 



Pc 



N 



(39) 



(40) 



(41) 



is a characteristic density. 

When performing the derivatives of the kinetic energy 
term, dimensional analysis along with the SVEA indicate 
that the second order derivative is smaller than the first 
by a factor 10^"*. We have therefore retained only the ad- 
vective contribution (first term on the right in Eq. (|36|) ) 
with a dimensionless velocity 



Lfll 2m 



(ki+k). 



(42) 



Finally, the rescaled atom-atom scattering strength is 



(43) 



The factor 7/(5l in Eqs. (|35|l and (|36() is due to the 
inclusion of a phenomenological decay term to account 
for losses due to spontaneous emission from the excited 
state. 

The evolution of the optical field is governed by the 
Maxwell wave equation coupled to the macroscopic po- 
larization of the condensate. Treating the pump laser as 
a continuous wave we find 
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d£L 



(1 - 1^1 5l) 



(44) 



■ 11 \^^^ 



(45) 



where ^ = 2:/^, C = ^/L and 



(46) 



Equations and (gSJ, together with Eqs. and 
(|35|l . fully describe the multiniode dynamics of the slowly 
varying envelopes of the electric and matter-wave fields. 
They may be solved analytically in the short-time regime 
as shown in Appendix A. In this regime the fields ^±k 
obey 



*fc(r,T) =^fe,(0)/o(2V7A) 



(47) 



where tpk' (0) is the square root of the recoiling mode 
density at r = and A is a function of ( defined in 
Appendix A. A similar result was previously obtained in 
Ref. JiJ. 



IV. NUMERICAL RESULTS 

A. General considerations 

This section presents selected results from numerical 
simulations of the onset and growth of superradiant scat- 
tering in a cigar-shaped condensate. We assume here that 
all atoms in the condensate at temperature T = are 
initially in the electronic ground state and approxi- 
mate their center-of-mass wave function by a separable 
Thomas- Fermi profile (in dimensionless units), 



i 

where the product is over the rescaled spatial coordinates 
TZi and Ci is the corresponding condensate width |26|. 
Close to the Thomas-Fermi radii, where the Thomas- 
Fermi approximation breaks down, we let the wavefunc- 
tion go to zero smoothly by matching both the function 
and its derivative to a Gaussian tail. This is necessary to 
prevent numerical instabilities resulting from the advec- 
tive term in Eq. (|36|l . To remain true to the experiments 
of Ref. (l| , we choose the aspect ratio of the condensate 
to be w/L Ki 0.1. 

As we have seen, the amplitudes of the end-fire modes 
at the onset of the classical regime are stochastic vari- 
ables the values of which must be selected at random 




FIG. 3: Typical time evolution of the EFM intensity. These 
simulations neglect dissipation, mean-field interactions and 
effects due to photon recoil. The curves correspond to: solid 
line N = 6.0, dashed line K = 3.0, gray line H = 1.0 and dot- 
dashed line N = 1.0 with 'y/Sh ~ 0.04. The x's mark points 
at which the spatial field and absorption profiles are plotted 
in Figs. 




FIG. 4: Density profile of the condensate at r = (top solid 
line) and of the atomic side mode "if+k at half the maximum 
EFM intensity, as indicated by x's in Fig. |3 solid line H = 
6.0, dashed line H = 3.0 and gray line K — 1.0. Mean-field 
interactions, dissipation and effects due to photon recoil are 
neglected. 



from run-to-run consistently with the results of Sec. II. 
However, to set the stage for our discussion, we consider 
first a simplified situation with only one left-recoiling and 
one right-recoiling side-mode of equal initial amplitudes 
small compared to the amplitude of ^'o, but with the 
same spatial structure. We include, for the moment, only 
the atom-light coupling and neglect the mean-field inter- 
action, photon recoil and dissipation. 

Figure |3| illustrates examples of superradiant EFM 
pulse shapes, for different values of the atom-field cou- 
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FIG. 5: Typical spatial envelopes of EFM electric fields at half 
maximum EFM intensity as indicated by x's in Fig. |21 solid 
line K — 6.0, dashed line H = 3.0, gray line H — 1.0. The peak 
in intensity for K = 1.0 and 3.0 arises due to scattering of EFM 
photons back into the probe beam. Mean-field interactions, 
dissipation and effects due to photon recoil are neglected here. 




FIG. 6: Absorption profile of the pump pulse along the z axis 
at half maximum EFM intensity as indicated by x's in Fig. |3] 
solid line K — 6.0, dashed line K = 3.0 and gray line K = 1.0. 
Mean-field interactions, dissipation and effects due to photon 
recoil are neglected here. 



pling strength H. To facilitate comparison we scaled 
the peak intensities to 1, and the time axes so that the 
peaks coincide. The gray line, corresponding to H = 1.0, 
demonstrates a typical main pulse followed by a sec- 
ondary peak, an effect known as ringing 0. Increasing 
the coupling strength to H = 3.0 raises the secondary 
peak maximum relative to the first peak, as shown by 
the dashed line. The black line, for which H = 6.0, rep- 
resents a qualitatively different regime: The superradi- 
ant maximum is now followed by a slowly decaying tail 
rather than ringing. Dissipation also leads to suppression 
of the ringing, as shown by the dot-dashed line which is 




FIG. 7: Time evolution of EFM intensity for an initial density 
of the recoiling modes 10~® (solid black line), 10~^ (dashed 
line) and 10~^ (gray line) times smaller than the condensate 
density. 



for j/Sl = 0.04 and H = 1.0. The appearance of ring- 
ing with increased coupling strength has been observed 
experimentally. Fig. 3(b) of Ref. 

Figures 0] and El show spatial profiles along the long 
axis z of the condensate, of the atomic side mode "if+k 
and of the corresponding EFM field £^k- Here the same 
values of N as in Fig. |3| were used: solid line H = 6.0, 
dashed line H = 3.0 and gray line H = 1.0. The EFM 
profiles have again been scaled for comparison. All z- 
axis profiles correspond to the times indicated by the x's 
in Fig. 131 In Fig. ^the atomic mode '^+k grows from 
the left edge of l^foKC)^, a consequence of the build-up of 
the EFM electric field as it moves across the condensate, 
see Fig. |S1 For clarity the mode '^-k is not shown. It 
is simply a mirror image of 'I'+fe. Interestingly, in the 
absence of dissipation the weaker coupling strength is 
more efficient at scattering atoms to ground state |2) over 
a single superradiant pulse, as evidenced by the gray line 
in Fig. ^ Keep in mind, however, that the time scale 
is slower for the weaker coupling, a feature not apparent 
due to the scaling of time in the figure. 

For the cases H = 1.0 and 3.0, the EFM profile peaks 
before reaching the edge of the condensate, see Fig. 
During the dynamical evolution this peak appears shortly 
after the first intensity maximum in Fig. (31 It is a con- 
sequense of the EFM field being scattered back into the 
probe beam from atoms in ground state |2). As pointed 
out in Ref. this may cause a strong EFM field to exist 
within the atomic system despite small emission outside 
the sample. 

The fraction of absorbed laser intensity 

is shown in Fig. j^lfor x ^ w after the laser field exits the 
condensate, again for the same values of H as in Fig. O 



8 




2 4 6 8 10 12 



FIG. 8: Comparison of EFM intensity with 7/5l = 8 x 10"^, 
H = 1.7 for superradiance only (black solid line), superra- 
diance and photon recoil (dashed line), superradiance and 
mean- field interaction (dot-dashed line), and superradiance 
with both mean-field interactions and photon recoil included 
(gray line). 



(a) 




0.5 


1 1.5 

t [ms] 


(b) A 





2 4 6 8 10 12 14 
T 



FIG. 9: Comparison of (a) experimental photomultiplier sig- 
nal of EFM intensity to, see 0), (b) a theoretical simulation 
with K = 1.7. Experimental data by D. Stamper-Kurn et al. 



and at half the EFM intensity maximum. The two cen- 
tral peaks in each absorption profile appear close to the 
edges of the condensate early in the time evolution and 
migrate towards each other. Secondary features, such as 
regions of gain around the edges of the condensate, ap- 
pear shortly after the EFM pulse maximum for H = 1.0 
and 3.0. This gain is a signature indicating that EFM 
light is being scattered back into the pump field. For 
H — 6.0 the growth of secondary absorption peaks at the 
edges of the condensate, as shown by the black line, is 
the dominant secondary feature rather than gain. 

The semi-classical model is valid provided the relevant 
modes are macroscopically occupied. In the above simu- 



lations the initial seed amplitude of the recoiling atomic 
modes was treated as a free parameter that may in gen- 
eral be chosen to fit experimental data. To study the 
effects of this choice of initial amplitude we plot in Fig.[7| 
the time evolution of EFM intensity for initial densi- 
ties of the recoiling modes being respectively 10"^ (solid 
black line), 10~^ (dashed line) and 10"'' (gray line) times 
smaller than the condensate density. For smaller initial 
amplitude the time to peak superradiance is delayed and 
the ringing peak is higher relative to the intensity maxi- 
mum. 

Now consider the effects of mean-field interactions and 
photon recoil. Figure |H1 shows the time evolution of the 
EFM intensity for 7/(5l = 8 x 10"^ ^nd H = 1.7. The 
results are for simulations taking into account superra- 
diance only (solid black line), superradiance and recoil 
velocity (dashed line), superradiance and mean-field in- 
teraction (dot-dashed line) , and superradiance with both 
mean- field interactions and recoil velocity (gray line). 
The recoil velocity was so chosen that an atom would tra- 
verse roughly half a condensate width over At = 12.0. 
The primary effect of including the photon recoil and 
mean-field interactions is to reduce the EFM amplitude. 
Note that the height of the ringing peak is increased 
somewhat relative to the first intensity maximum when 
including only mean-field interactions. But generally 
speaking we do not find significant qualitative changes 
in the spatial and temporal evolution of either matter or 
light fields other than a reduction in overall amplitudes 
and slight relative delays. 



B. Comparison with experiment 

The primary purpose of the multimode treatment is 
to describe the quantum fluctuations that govern the ini- 
tial evolution of the system. As discussed in Sec. Ullthis 
entails randomly choosing seed amplitudes for the intial 
state of the recoiling atomic modes with the appropriate 
statistical properties. We have found, however, that the 
multimode description gives the same results as a single- 
mode model as far as EFM intensity and absorption pro- 
files are concerned, provided that the sum of initial seed 
densities in the multimode case is equal to the seed den- 
sity of the single-mode system. For a single run of the 
experiment a single-mode description is therefore suffi- 
cient if the seed densities are chosen judiciously. We thus 
turn to a comparison of the single-mode theory to data 
obtained in experiments at UC Berkeley 

As described in Ref. 1], a Bose condensate of = 
1.6 X 10^ ®^Rb atoms trapped in an loffe-Pritchard trap 
with trap frequencies of uj^.y.z = 27r(48, 48, 5)s~^ was su- 
perradiantly pumped via the _F = 1^_F' = 1D1 transi- 
tion from the |F = ^,mf = —1) to the \F = 2,m-f = 1) 
groundstates via a sequence often pump pulses of 100 fis, 
each separated by a 68 fj,s delay. Fig. Ela) shows a typical 
EFM pulse sequence recorded on a photomultiplier tube 
(PMT) . The rise time between pulses is a consequence of 
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FIG. 10: (a) Experimental coherent absorption imaging data, see in which darkest regions indicate strongest absorption, 
with a maximum absorption of ~ 15%. (b) Simulated absorption images for — 1.7 and an initial density of recoiling modes 
IvI'oP/lO^. (c) Comparison of experimental (circles) to simulated (gray lines) transverse averaged absorption profiles. The 
maximum in the experimental data is ~ 6%. In all figures time evolves from top to bottom and the horizontal frame width is 
L ~ 200fim. Experimental data by D. Stamper-Kurn et al. [24|. 




FIG. 11: Example of (a) Typical absorption image, see Q, 
and (b) its averaged profile, demonstrating asymmetry that 
may arise from quantum fluctuations. Experimental data by 
D. Stamper-Kurn et al. [2^ . 



the slow response time of the PMT. The lack of ringing 
indicates that dissipation plays a significant role in the 
dynamics. Dissipation due to spontaneous emission may 
be estimated using ^/5l = oj'^cP /Sneohc^dL « 8 x lO^'^ 
|25| . Given an s-wave scattering length for *^Rb of 
a « IO^Oq where Qq is the Bohr radius, we find H ~ 1.7. 
The experimental data agree qualitatively with the three 
dimensional simulation, see Fig. Etb). Here the initial 
density of recoiling atoms was chosen 10"* times smaller 
than the condensate density, and the dimensionless time 
scale was such that an atom would traverse roughly half 
a condensate width over the course of the pulse sequence 
as in the Berkeley experiments. The PMT response was 
accounted for by convoluting the EFM intensity with a 
response function with response time of the order of the 
pulse separation. 



Experimental absorption images corresponding to the 
PMT signal of Fig. IHIa) are shown in Fig. Iior al. where 
time evolves from top to bottom and each image is the 
time averaged absorption over the duration of each pulse 
in the sequence. Consistently with the observation that 
in this experimental run the asymmetry in the absorption 
profile is small, we have taken the initial amplitudes of 
the atomic side modes 'i>±k to be equal. The resulting ab- 
sorption images are shown in Fig. llOf b). An effect not 
resolved in Ref. is the weak absorption seen at the 
upper- and lower edges of the condensate in the lower 
six images of Fig. Iiur bl. This absorption arises since 
the edges of the condensate are less dense than the cen- 
ter, thus delaying superradiance in time and intensity. 
A comparison of our simulation to transverse averaged 
experimental absorption profiles is shown in Fig. Iiur c). 
The qualitative agreement is quite satisfactory although 
the value of the peak absorption is underestimated by 
a factor of '--^ 2. Note that the averaging procedure su- 
presses gain and other secondary features seen in the sim- 
ulation of Fig. [TUT b). thus leading to the simple double 
peaked structure of Fig. Iiur c). The simulation further- 
more suggests that these features are too small in spatial 
extent and amplitude to be resolved in the experimen- 
tal data. We emphasize that the only free parameters in 
this simulation are the initial seed amplitude and the to- 
tal time duration of the pulse sequence. They are tightly 
constrained by, respectively, the rise time to peak super- 
radiance relative to the EFM pulse width and the EFM 
pulse width relative to the overall time duration of the 
pulse sequence. 
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FIG. 12: Absorption profile A{z) of the input laser field. Figures (a)-(b) are representative profiles resulting from single 
realizations of a superradiance experiment for geometries with (a) m = 5 and (b) m = 20 quasi-modes each in the left- and 
right EFM directions. The peak separation, Aa, as indicated by the dotted lines in (a) characterizes the asymmetry. Figures 
(c) and (d) give the frequency of occurrence of the absorption profiles resulting from 500 realizations for (c) m — 5 and (d) 
m = 20. As indicated by the gray scale in the color bar, the darkest shade of gray corresponds to the most frequent occurrences, 
(e) Distribution of the asymmetry of absorption peaks from 500 realizations, for systems with m = 5 (black line), m = 10 (gray 
line) and m = 20 (dashed line) quasi-modes in both the left- and right EFM directions. 



C. Shot-to-shot fluctuations 

We now turn to the shot-to-shot fluctuations result- 
ing from the quantum noise-dominated early stages of 
the emission process. In Fig. Illf a) we show a typical 
absorption image obtained experimentally and exhibit- 
ing strong asymmetry, as well as its transverse averaged 
proflle in Fig. Illf b). To account for this asymmetry 
we solve the multimode Eqs. (jSnj-lISnil and |g3I)-gSJl for 
H = 3.0, neglecting mean-field interactions and effects 
due to photon recoil, with the initial amplitude of each 
quasi-mode chosen from the probability distribution 
independently of the other modes. The average particle 
number fiq is taken to be the same for all distributions. 

Figures Efa)-(b) show absorption profiles for two ran- 
domly chosen realizations of a superradiance experiment 
for systems with (a) m — 5 and (b) m = 20 quasi-modes 
each in the left- and right EFM directions of propagation. 
In all cases the absorption profile is shown at the time 
of the first EFM intensity maximum in either the left- 
or right traveling EFM. The separation AA, as indicated 
by the dotted fines in (a), may be used to characterize 
the asymmetry. Figure [T^ c)-fd) shows a gray scale plot 
constructed from 500 realizations of the superradiance 
experiment for (c) m — 5 and (d) m = 20 quasi-modes 
each in the left- and right travelling EFMs. Each pro- 
file is again taken at the first EFM intensity maximum. 
The gray scale represents the frequency of occurrence of 
absorption, with darkest shade of gray being the highest 
frequency. 

The distribution of absorption asymmetries, AA, for a 



random set of 500 realizations of the experiment is plot- 
ted in Fig. I12r e) for m = 5 (black line), m = 10 (gray 
line) and m = 20 (dashed line) quasi-modes each in the 
left- and right EFM directions. As may be expected the 
width of the distribution decreases with increasing num- 
ber of modes. 



V. CONCLUSION 

We have presented a multimode model of Raman su- 
perradiance scattering to explain recent experiments that 
image Bose-condensates based on the abovementioned 
principle. The superradiant scattering into end-fire- 
modes is seen to lead to non-trivial, time-dependent, spa- 
tial structure in both matter- and photon fields. The 
microscopic quantum fluctuations that trigger the initial 
superradiant dynamics lead to macroscopic fluctuations 
and asymmetries in the spatial features of absorption im- 
ages during later stages. The modeled absorption profiles 
are in good qualitative agreement with experimental ob- 
servations, confirming that superradiant scattering can 
be used to probe the condensed phase of a BEC at fi- 
nite temperature while remaining insensitive to the non- 
condensed phase. 
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APPENDIX A 

The short-time expressions may be obtained by making 
an approximation of zero depletion for the laser field, 
fi, ~ 1, and keeping only highest order terms in the 
equations of motion for the atomic fields and end-fire- 
modes, 

I^Lp^o (A-1) 

^k^i^-o (A-2) 

*LfefL*o (A-3) 

Equation ljA-l|) may immediately be solved to give 

^'o(t) = ■^'^^e'^^e-"^ (A-4) 

where rj = 



d*o(r,i) 




i 


dT 
d*k(r,i) 


(1 


- il/S) 
i 


dT 

sign(fc)— ? 


^ (1 
(1 


- il/S) 



(T+(W) ' " ^""^ ^0*" 

condensate wave function which is assumed to be of the 
form ^Q°^ = (/)^(/)(^ X)- Here T = y/L is the dimensionless 
length in the y-direction. We also assume that in the 
initial state the recoiling atomic modes have the same 
wavefunction as the condensate, but smaller in amplitude 
by a factor /3 < 1. At time t = the EFM field is, by 
direct integration of Eq. ljA-3p , 

sign(fc)ffc(C) - ^{iii~v)P\<^^^^r)?(^j\\<Pc?dC' + 



H(«/-i.)/3|0(5^x)|'A(C)/a, 



(A-5) 



where £q is chosen so that the boundary condition 
£_fe(l/2) = is satisfied and 



A(C) = a 



j\ \^c?dC + £!)j 



where 



(A- 



(A-7) 



Now substitute Eq. (|A-4|) into Eq. IjA-Sp and define 
= E'if^e"''^^'^'^ . Then upon taking the time derivative 
of Eq. I)A-3|1 and substitution of Eq. \A-2\ therein we 
find 



sign(fc) 



dTdQ 



(A- 



To solve this differential equation we first treat the case of 
zero dissipation, i.e. 77 — > 1 and u ^ using the method 



of Laplace Transforms. Taking the Laplace transform 
with respect to time of Eq. HA-8|I we obtain 



s — zNp Wq ' 



dC 



(A-9) 



in which E*_f.{(,.s) = C {^-kj- Equation (|A-9p is a lin- 
ear first order differential equation and my be solved by 
standard means. We find 

E*_,(C, s) = be^'^O/s _ ,H/3|0(^^x)lV«- (A-10) 

The EFM field is found by inversion of the Laplace trans- 
form 



1 p^i+ica 



(A-11) 



fi— too 



The contour of integration is chosen so as to include all 
poles of the integrand and the integral may be evaluated 
using the residue theorem. This is done by expanding 
both exponentials in the integrand in series and multiply- 
ing out the result. Integrating term- by-term the residue 
is found to be of the form 



where Iq is a modified Bessel function of the first kind of 



^ (n- l)\nl 

n=l ^ ' 

is a I 
order q. Thus 

£-l,(C,r) =1/3^^/1 (2V7a) 
and by substitution into Eq. l|A-3p 



(A-12) 



i/=0 



i/=0 



(A-13) 



(A-14) 



We now turn to the case with non-zero dissipation. 
Since the solution to the latter must give rise to Eq. 
(|A-13|I for — > 0, we suggest a solution of the form 



flfc ai(T)A -f a2(T)A2 ^ az{T)l\^ ^ .. 



(A-15) 



were the a„ (r) are time dependent coefficients that must 



vanish at r = and reduce to 



for 



0. Substi- 



tuting Eq. ljA-15|) into Eq. IjA-Sp one obtains a recursion 
relation for the nth coefficient 



da-a 



To satisfy the — > condition wc find that 



1 



(1 



-2vt\ 



2"-ij/"-i(n - l)!n! 
The latter leads to the desired solutions 

2z.A(C) 



(A-16) 



(A-17) 



r,(r,r)=z/3(l-z7/J) 



1 



h (giO) , 

(A-18) 



12 



and 

*,(r,T)=/3*(°)7o(5(C))- (A-19) 
The argument of the Bessel functions is given by 

P(0=(^^S)(1^)*. ,A..o) 

Finally, to find absorption profiles of the laser field we 
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/(C) = H (i77 - v) cP^ j " 4>{i,y)di. (A-23) 



For a Thomas-Fermi profile (pQ — (1/2)^ — C"^ and 

A(C) = cx (\{l/2f - {l/2fC + \e) . (A-24) 



